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Abstract 

We provide a solution to the problem of determining whether a target pure state can 
be asymptotically prepared using dissipative Markovian dynamics under fixed locality 
constraints. Beside recovering existing results for a large class of physically relevant 
entangled states, our approach has the advantage of providing an explicit stabilization 
test solely based on the input state and constraints of the problem. Connections with 
the formalism of frustration-free parent Hamiltonians are discussed, as well as control 
implementations in terms of a switching output-feedback law. 

1 Introduction 

While uncontrolled couplings between a quantum system of interest and its surrounding 
environment are responsible for unwanted non-unitary evolution and decoherence, it 
has also been long acknowledged that suitably engineering and exploiting the action 
of the environment may prove beneficial in a number of applications across quantum 
control and quantum information processing (Poyatos et al. 1996, Beige et al. 2000, 
Lloyd & Viola 2001). It is well known, in particular, that open-system dynamics 
are instrumental in control tasks such as robust quantum state preparation and rapid 
purification, and both open-loop and quantum feedback methods have been extensively 
investigated in this context (Combes et al. 2008, Wiseman & Milburn 2009, Ticozzi & 
Viola 2009, Schirmer & Wang 2010, Combes et al. 2010), including recent extensions 
to engineered quantum memories (Pastawski et al. 2010) and pointer states in the 
non-Markovian regime (Khodjasteh et al. 2011). 

Remarkably, it has also been recently shown that it is in principle possible to de- 
sign dissipative Markovian dynamics so that non-trivial strongly correlated quantum 
phases of matter are prepared in the steady state (Diehl et al. 2008) or the output of 
a desired quantum algorithm is retrieved as the asymptotic equilibrium (Verstraete et 
al. 2009). From a practical standpoint, scalability of such protocols for multipartite 
systems of increasing size is a key issue, as experimental constraints on the available 
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control operations may in fact limit the set of attainable states. Promising results have 
been obtained by Kraus et al. (2008) for a large class of entangled pure states, showing 
that Markovian dissipation acting non-trivially only on a finite maximum number of 
subsystems is, under generic conditions, sufficient to generate the desired state as the 
unique ground state of the resulting evolution. As proof-of-principle methodologies 
for engineering dissipation are becoming an experimental reality (Barreiro et al. 2011, 
Krauter et al. 2011), it is important to obtain a more complete theoretical characteri- 
zation of the set of attainable states under constrained control resources, as well as to 
explore schemes for synthesizing the required dissipative evolution. 

Building on our previous analysis (Ticozzi & Viola 2008, 2009, Ticozzi et al. 2010), 
in this work we address the problem of determining whether a target pure state of a 
finite-dimensional quantum system can be prepared employing "quasi-local" dissipative 
resources with respect to a fixed locality notion (see also Yamamoto's contribution to 
this volume for related results on infinite-dimensional Markovian Gaussian dissipation) . 
We provide a stabilizability analysis under locality-constrained Markovian control, in- 
cluding a direct test to verify whether a desired entangled pure state can be asymp- 
totically prepared. In addition to recovering existing results within a system-theoretic 
framework, our approach has the important advantage of using only two inputs: the 
desired state (control task) and a specified locality notion (control constraints), with- 
out requiring a representations of the state in the stabilizer, graph, or matrix-product 
formalisms. 

2 Problem definition and preliminary results 
2.1 Multipartite systems and locality of QDS generators 

We focus on quantum evolutions driven by a (time-independent) Markovian Master 
Equation (MME) (Gorini et al. 1976, Lindblad 1976, Alicki & Lendi 1987) in Lindblad 
form (ft = 1): 

p{t) = C{p{t)) = -i[H, p{t)\ + J2 {LkP{t)Ll ~ \{LlLk, p{t)}) , (1) 

k 

specified in terms of the Hamiltonian H = H'' and a finite set of noise (or Lindblad) 
operators {Lk}. We are interested on the asymptotic behavior of MMEs in which 
the operators H,{Lk} satisfy locality constraints. More precisely, let us consider a 
multipartite system Q, composed of n (distinguishable) subsystems, labeled with index 
a = 1, . . . , n, with associated da-dimensional Hilbert spaces Ha- Thus, T-Lq = ®a=i Ha- 
Let QS('H) and D{'H) denote the sets of linear operators and density operators on "H, 
respectively. It is easy to show (see e.g. Ticozzi & Viola 2008, proof of Theorem 2) 
that the semigroup generated by Eq. ([1} is factorized with respect to the multipartite 
structure, that is, the dynamical propagator 

n 

7^ -e^* = (g)ra,t Vt>0, 

with Ta,t Si CPTP map on %{'Ha), if and only if the following conditions hold: 

(i) Each Lk acts as the identity on all subsystems except (at most) one; 

(ii) H — Ha, where each Ha acts as the identity on all subsystems except (at 
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most) one. 

This motivates the following definitions: we say that a noise operator Lk is local if it 
acts as the identity on all subsystems except (at most) one, and that a Hamiltonian 
H is local if it can be written as a sum of terms with the same property. However, 
it is easy to verify that if a semigroup associated to local operators admits a unique 
stationary pure state, the latter must necessarily be a product state. Thus, in order 
for the MME ([TJ to admit stationary entangled states, it is necessary to weaken the 
locality constraints. 

We shall allow the semigroup dynamics to act in a non-local way only on certain 
subsets of subsystems, which we call neighborhoods. These can be generally specified 
as subsets of the set of indexes labeling the subsystems: 

Nj<Z{l,...,n}, j = l,...,M. 

In analogy with the strictly local case, we say that a noise operator L is Quasi-Local 
(QL) if there exists a neighborhood Nj such that: 

where Lj\f- accounts for the action of L on the subsystems included in A/}, and /j^-^. := 
^ai^M identity on the remaining subsystems. Similarly, a Hamiltonian is QL 

if it admits a decomposition into a sum of QL terms: 

3 

A MME will be called QL if both its Hamiltonian and noise operators are QL. It is 
well known that the decomposition into Hamiltonian and dissipative part of ([1} is not 
unique: nevertheless, the QL property remains well defined since the freedom in the 
representation does not affect the tensor structure of H and {Lk}- The above way of 
introducing locality constraints is very general and encompasses a number of specific 
notions that have been used in the physical literature, notably in situations where the 
neighborhoods are associated with sets of nearest neighbors sites on a graph or lattice, 
and/or one is forced to consider Hamiltonian and noise generators with a weight no 
larger than t (so-called t-body interactions), see also Kraus et al. 2008, Verstraete et 
at 2009. 

We are interested in states that can be prepared (or, more precisely, stabilized) by 
means of MME dynamics with QL operators. Recall that an invariant state p for a 
system driven by ([!]) is said to be Globally Asymptotically Stable (GAS) if for every 
initial condition po we have 

lim e^*(po) = P- 

t—>oo 

In particular, following Kraus et al. (2008), the aim of this paper is to characterize 
pure states that can be rendered GAS by purely dissipative dynamics, for which the 
state is "dark". More precisely: 

Definition 1 A pure state pd = G 'S{'Hq), is Dissipatively Quast-Locally Sta- 

bihzable (DQLS) if there exist QL operators {Dk}k^i,...,K on Hq, with D^l'i') ~ 0, for 
all k and Dk acting non-trivially on (at most) one neighborhood, such that pd is GAS 
for 

P-^n[p] = Yl {^^P^l - \{DlD>^' P})- (2) 
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We will provide a test for determining whether a state is DQLS, and in doing so, we 
will also show how assuming a single QL noise operator for each neighborhood does 
not restrict the class of stabilizable states. From now on, we thus let K = M, and 

We begin by noting that if a pure state is factorized, then we can realize its tensor 
components "locally" with respect to its subsystems (see Ticozzi & Viola 2008, Ticozzi 
et al. 2010) for stabilization of arbitrary quantum states in a given system with simple 
generators, involving a single noise term). Thus, we can iteratively reduce the problem 
to subproblems on disjoint subsets of subsystems, until the states to be stabilized are 
either entangled, or completely factorized. A preliminary result is that the DQLS 
property is preserved by arbitrary Local Unitary (LU) transformations, of the form 
U — ^2=1 -'-'1 order to show this, the following lemma is needed. 

Lemma 1 Let C denote the Lindblad generator associated to operators H,{Lk}- Then 
for every unitary operator U we have 

U£.[U''pU]U'' = C'[p] (3) 

where C' is the semigroup generator associated to H' = UHU\ L'^. = ULkU\ and 

l7e'^*[(7Vof/]t^^ = e^'*M, > 0. (4) 

Identity ([2} is easily proven by direct computation, while Q follows directly from the 
properties of the (matrix) exponential. The desired invariance of the QL stabilizable 
set under LU transformation follows: 

Proposition 1 If p is DQLS and U is LU, then p = U pU'^ is also DQLS. 

Proof. Assume that the generator C associated to QL operators {Dk} stabilizes p. Since 
p is GAS, for any initial condition po we may write 

p' = lim Ue'^'[po]U'' = lim U{e^'[U'' poU])U\ (5) 

t— >-|-00 t—¥-{-00 

By applying Lemma [Jl it suffices to show that each Dj. — UDkU'' is QL. Since Dk = 
Dm^ ® Im^, we have D'^ = U{Dm^ ® ^^rJ'^'' = {Um^Dm^UI;^^ ® I^^, where Um^ ■= 
^leAfk Hence D'^. is QL. □ 



2.2 QDS for unconstrained stabilization 

We next collect some stabilization results that do not directly incorporate any locality 
constraint, but will prove instrumental to our aim. Let Hs :~ span{|'I')}. Given 
Corollary 1 in Ticozzi & Viola 2008, pd is invariant if and only if 

k 

where we have used the natural block representation induced by the partition Hq — 
Hs ® Hs Ei'iid labeled the blocks as 



■ Xs 


Xp ' 




Xr 



Assume pd to be invariant. Hence must be a common eigenvector of each Lk. 
Call the corresponding eigenvalue lk = Ls,k. By Lemma 2 in Ticozzi & Viola (2008), 



Lk 



^S.k 



4 



the MME generator is invariant upon substituting Lt with Lk — Lk — Ikl, and H = 
H + i^k^Pk^k — f-kL\.). In this way, we have Ls,k = for all fc, so that Hp must be 
zero in order to fulfill the above condition. Thus, H is block-diagonal, with \^) being 
an eigenvector with eigenvalue h = Hs- Using this representation for the generator, 
we can let Ls.k = tk ~ 0, and Hp = = Hq. This further motivates the use of noise 
operators Dk such that Dk\'i') = in the DQLS definition. 

Lemma 2 An invariant pd = |^)(^| is GAS for the MME ^ if there are no invariant 
common (proper) subspaces for {Lk} other that Hs ~ span{|^')}. 

Proof.By Lemma 8 and Theorem 9 in Ticozzi & Viola 2009, pd is GAS if and only if 
there are no other invariant subspaces for the dynamics. Given the conditions on Lk for 
the invariance of a subspace, pd is GAS as long as there are no other invariant common 
subspaces for the matrices Lk- □ 

Based on the above characterization, in order to ensure the DQLS property it 
suffices to find operators {Dk}k=i C such that Hs is the unique common 

(proper) invariant subspace for the {Dk}- 



3 Characterization of DQLS states 
3.1 Main result 

Our main tool for investigation will be provided by the reduced states that the target 
state Pd induces with respect to the given local structure. Let us define: 

pJ^k = traceyg-Jpd), (6) 

where tracej^^ indicates the partial trace over the tensor complement of the neighbor- 
hood A/fc, namely Hj^^ = ®a<^j^,, '^'^'^ following Lemma follows from the properties 
of the partial trace: 

Lemma 3 supp{pd) C f]^ supp{pj^f^ (g) Ij^J. 

Proo/.From the spectral decomposition pAA*, ~ X^gP^n^, we can construct a resolution 
of the identity {II,} such that pAA^ ® I//^ = SgP^n, (g) I//^, where, by definition 
of the partial trace, pq = tTa.ce{pdTlq ® I//,,)- Pq =0 for some q, then it must 
be Pd(n^ ig) /jv-j.) = and therefore supp(pd) -L supp(n^ g) -Ta/-^). Thus, supp(pti) C 
[jq supp{pqUq (» I^fJ = supp(paaJ, for all k. □ 

Let us now focus on QL noise operators Dk ~ Dj^^ ® Ifj-^ such that Dk\^) ~ 0. 

Lemma 4 Assume that a set {Dk} makes pd ~ DQLS- Then, for each k, we 

have supp{pj\{^ ) C ker(DArjj ) . 

Proof-Since is by hypothesis in the kernel of each Dk, with respect to the decom- 
position Hq = Hs ffi Hs every Dk must be of block form: 



Dk 



Dp.k 
Duk 



which immediately implies DkPdDj, = 0. It then follows that tracej^-^ (-Dfep^-D^) = = 
trace^^(DAAfc <8) Ij^^pdD]^^ ® ^aTj,)- Therefore, it also follows that D^^pM^D'^j^^ = 0. 
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If we consider the spectral decompositon pAf^ = X]j 1j\4'j) {4'j\y with qj > 0, the latter 
condition implies that, for each j, Dj^^\(j)j){(f)j\D]^^ — 0. Thus, it must be supp(pA/'fc) CI 
ker(E)Ark), as stated. □ 

Theorem 1 A pure state pd ~ \tp){tp\ is DQLS if and only if 

supp{pd) = Pi supp{pM^ ® /yg-J = Pi Hjv-fc . (7) 

k fe 

Proo/.Given Lemmas 3 and 4, for any set {Du} that make pd DQLS we have: 

supp(pd) C p supp(pAAfe 8) /atJ ^ P ker(DAAk ® IatJ- 
fe fe 

By negation, assume that supp(pd) $ Dfe suPP(PA/'fc (8) Im^)- Then there would be (at 
least) another invariant state in the intersection of the kernels of the noise operators, 
contradicting the fact that pd is DQLS. Thus, a necessary condition for pd to be GAS 
is that supp(pti) = plj. supp(pA/'fc ® ^Mk^- Conversely, if the latter condition is satisfied, 
then for each k we can construct operators D^f^, that render each supp(pArj,) GAS on 
"HATfc (see e.g. Ticozzi & Viola 2008, 2009, Ticozzi et al. 2010 for explicit construc- 
tions). Then ker(DA/ij ®'^fj'^) = supp(pA4), and there cannot be any other invariant 
subspace. By Lemma [21 Pd is hence rendered GAS by QL noise operators. □ 



3.2 An equivalent characterization: QL parent Hamiltoni- 
ans 

Consider a QL Hamiltonian H — "^^k Hk, Hk = H^f|^ '^^JVk ■ ^ Pure state pd = 
is called a frustration-free ground state if 

= min A(//fe), Vfc, 

where A(-) denotes the spectrum of a matrix. A QL Hamiltonian is called a parent 
Hamiltonian if it admits a unique frustration-free ground state (Perez-Garcia et al. 
2007). 

Suppose that a pure state admits a QL parent Hamiltonian H. Then the QL 
structure of the latter can be naturally used to derive a stabilizing semigroup: it suffices 
to implement QL operators Lk that stabilize the eigenspace associated to the minimum 
eigenvalue of each Hk- In view of Theorem 1, it is easy to show that this condition is 
also necessary: 

Corollary 1 A state l^*) is DQLS if and only if it is the ground state of a QL parent 
Hamiltonian. 

Proo/. Without loss of generality we can consider QL Hamiltonians H = ^f. Hk, where 
each Hk is a projection. Let pd be DQLS, and define H^ := Hj^^ ® /jv'^i with 11;^^ 
being the orthogonal projector onto the orthogonal of the support of PaAj,, that is, 
T-Lj4^Qsupp{pj^^) . Given Theorem[T] l^*) is the unique pure state in supp(pA/'fc®Jjv'fc)i 
and thus the unique state in the kernel of all the H^. Conversely, if a QL parent 
Hamiltonian exists, to each Hk we can associate an Lfe that asymptotically stabilizes 
its kernel. A single operator per neighborhood is in principle always sufficient (see 
Ticozzi & Viola 2008, 2009, Ticozzi et al. 2010 for explicit constructions and examples 
of Lk stabilizing a desired subspace). □ 
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The above result directly relates our approach to the one pursued in Kraus et al, 
2008, Verstraete et al. 2009, and a few remarks are in order. In these works it has 
been shown that Matrix Product States (MPS) are QL stabilizable, up to a condition 
(so-called injectivity) that is believed to be generic (Perez-Garcia et al. 2007). MPS 
states that allow for a compact representation (that is, in the corresponding "valence- 
bond picture", those with a small bond dimension) are of key interest in condensed 
matter as well as quantum information processing (Verstraete al. 2006, Perez-Garcia 
et al. 2007, Perez-Garcia et al. 2008). However, any pure state admits a (canoni- 
cal) MPS representation if sufficiently large bond dimensions are allowed, suggesting 
that arbitrary pure states would be DQLS. The problem with this reasoning is that 
the locality notion that is needed in order to allow stabilization of a certain MPS is 
in general induced by the state itself. The number of elements to be included in each 
neighborhood is finite but need not be small: while this is both adequate and sufficient 
for addressing many relevant questions in many-body physics (where typically a ther- 
modynamically large number of subsystems is considered), engineering the dissipative 
process may entail interactions that are not easily available in experimental settings. 
For this reason, our approach may be more suitable for control-oriented applications. 
It is also worth noting that the injectivity property is sufficient but not necessary for 
the state to admit a QL parent Hamiltonian (an example on a two-dimensional lattice 
is provided in Pcrcz-Garcia et al. 2007). Once the locality notion is fixed, our test for 
DQLS can be performed irrespective of the details of the MPS representation, and it is 
thus not aifected by whether the latter is injective or not (rather, our DQLS test may 
be used to output a QL parent Hamiltonian if so desired). 

3.3 Examples 

• GHZ-states and W-states.— Consider an n-qubit system and a target GHZ state 
PGHZ = |*)(*|, with 1*) = |4'ghz) = (I000...0) |lll...l))/\/2. Any reduced 
state on any (nontrivial) neighborhood is an equiprobable mixture of |000 ... 0) and 
I 111 ... 1). It is then immediate to sec that 

span{|000 ... 0), |111 . . . 1)} C f|supp(pAr, ® 

k 

and hence pghz is not DQLS. In a similar way, for any n the W state pw = |*)(*|, 
with I*) = |*w) = (|100 ...0) + |010 ... 0) -h ... -h |000 . . . l})/^/n has reduced states 
that are statistical mixtures of |000 ... 0) and a smaller W state |^w')i of the dimension 
of the neighborhood. Thus, 

span{|000 ... 0), |*w')} ^ f] supp{pM^ (g) 7^J, 

k 

and Pw is not DQLS (except in trivial limits, see also below). Note that for arbitrary n, 
both PGHZ and pw are known to be (non-injective) MPS with (optimal) bond dimension 
equal to two. 

• Stabilizer and graph states.— A large class of states does admit a QL description, 
and in turn they arc DQLS. Among these are stabilizer states, and general graph 
states. Here the relevant neighborhoods are those that include all the nodes connected 
to a given one by an edge of the graph. The details are worked out in Kraus et al. 
2008. Notice that GHZ states are indeed graph states, but only associated to star (or 
completely connected) graphs. Hence, relative to the locality notion naturally induced 
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by the graph, any central node has a neighborhood which encompasses the whole graph, 
rendering the constraints trivial. 

• DQLS states beyond graph states.— Consider a 4-qubit system arranged on 
a linear graph, with (up to) 3-body interactions. The two neighborhoods A/i — 
{l,2,3},A/'2 = {2,3,4} are sufficient to cover all the subsystems, and contain all the 
smaller ones. The state pr = with 

I*) = |*t) = (|1100) + 11010) + 11001) + 10110) + 10101) + 10011))/V6, 

is not a graph state, since if we measure any qubit in the standard basis, we are left with 
W states on the remaining subsystems, which are known not to be graph states. In 
contrast. Proposition 9 of Hein et al. 2004 ensures that the conditional reduced states 
for a graph state would have to be graph states as well. Nonetheless, by constructing 
the reduced states and intersecting their supports one can establish directly that |^')t 
is indeed DQLS. 



4 Switched feedback implementation 



From Theorem[T]it follows that a DQLS state can be asymptotically prepared provided 
we can engineer QL noise operators Dk — Dj\f^ (g) Ij^ that stabilize the support of 



with the blocks Dp,k, Dn^^ such that the support of paA^ is 



each reduced state p^,. on each neighborhood. Restricting to Hat^ , we must have 
Dp,k 
DR,k 

attractive, that is, such that no invariant subspace is contained in its complement. 
Following the ideas of Ticozzi & Viola (2009), Ticozzi et al. (2010), a natural explicit 
choice is to consider noise operators with the following structure: 



Dp.k = 





: 

h 



■0 £2 

£3 ■•• 



(8) 



If the above QL Lindblad operators are not directly available for open-loop im- 
plementation, a well studied strategy for synthesizing attractive Markovian dynamics 
is provided by continuous measurements and output feedback. In the absence of ad- 
ditional dissipative channels, and assuming perfect detection, the relevant Feedback 
Master Equation takes the form (Wiseman & Milburn 2009): 

m ^-i[H + H,+ \{FM + M^F), p(t)] + Lfp{t)L\ - \ {Lt.L/, p(t)} , 

where He is a time-independent control Hamiltonian, F = F^ and M denote respec- 
tively the feedback Hamiltonian and the measurement operator, and Lf := M — iF. 
Necessary and sufficient conditions for the existence of open- and closed- loop Hamil- 
tonian control that stabilizes a desired subspace have been provided in Ticozzi & Viola 
2008, 2009. 

In order to exploit the existing techniques in the current multipartite setting, it 
would be necessary to implement measurements and feedback in each neighborhood. 
If the measurement operators do not commute, however, one would have to carefully 
scrutinize the validity of the model and the consequences of "conflicting" stochastic 
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back-actions when acting simultaneously on overlapping neighborhoods. These difficul- 
ties can be bypassed by resorting to a cyclic switching of the control laws. Consider a 
DQLS state pd and the family of generators {£k}kLi, Cklp] = D^pD], — ^{D\Dk, p}, 
with Dk such that supp(pA^j. (g) Iff,,) is the unique invariant subspace for Ck- Define a 
switching interval r > and the cyclic switching law j{t) — [t/rM] + 1. We can then 
establish the following: 

Theorem 2 There exists QL {Dt} such that pd is GAS for the switched evolution 

Proof. Consider the trace-preserving, completely-positive maps Tj{p) = e^^'^[p]. It is 
easy to see that pd is invariant for each 7j : as a corollary of Theorem 1 in Bolognani 
& Ticozzi 2010, it follows that pd is GAS if it is the only invariant state for T — 
Tm o • ■ ■ o 71 . Assume that p is invariant for T: then either it is fixed for all Tk , which 
means that necessarily p = pd, or there exists a periodic cycle. Since each Tj is a trace- 
distance contraction (Alicki & Lendi 1987), this means that each map preserves the 
trace distance, that is, \\Tj{pd — p)\\i ~ \\pd — p\\i. This would in turn imply that each 
Tj admits eigenvalues on the unit circle, and hence each jCk would have imaginary ones. 
However, if we choose Dk as in Eq. ((S)), in vectorized form the Liouvillian generator 
reads 

jCk = £>r ® £>fc - ® DlDk - ^{DlDkf ® /, 

which is an upper triangular matrix with eigenvalues either equal to zero or { — (^j + 
£f)/2}. Therefore, for this choice pd is the only invariant pure state state for T and 
hence it is GAS. □ 



5 Concluding remarks 

We have presented a characterization of DQLS pure states for fixed locality constraints, 
from a control perspective. As a byproduct of our main result, an easily automated 
algorithm for checking DQLS states is readily devised. The necessary steps entail: (1) 
calculating the reduced states on all the neighborhoods specifying the QL notion; (2) 
computing their tensor products with the identity on the remaining subsystems, and 
the relative supports; (3) finding the intersection of these subspaces. If such intersec- 
tion coincides with the support of the target state alone, the latter is DQLS. If so, we 
have additionally showed that the required Markovian dynamics can in principle be im- 
plemented by switching output-feedback control. While we considered homodyne-type 
continuous-time feedback MME, the study of discrete-time strategies is also possible 
along similar lines, see also Bolognani & Ticozzi (2010), Barreiro et al. (2011). 

Our present results have been derived under two main assumptions: the absence 
of underlying free dynamics, and the use of purely dissipative control (no Hamiltonian 
control involved). In case a drift internal dynamics is present, the same approach 
can be adapted to determine what can be attained by dissipative control. When we 
additionally allow for Hamiltonian control, one may employ the algorithm described 
in Section III.B of Ticozzi et al. (2011) to search for a viable QL Hamiltonian when 
dissipation alone fails. Nonetheless, in the presence of locality constraints a more 
efficient design strategy may be available: an in-depth analysis of these issues will be 
presented elsewhere. 

It is also worth noting that in various experimental situations the available dissi- 
pative state preparation procedures involve two steps: first, enact local noise operators 
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that prepare a known pure state that is factorizcd; next, use open-loop coherent control 
to steer the system on the desired entangled target. The approach we discussed here is 
believed to have an advantage in terms of the overall robustness against initialization 
errors and finite-time perturbations of the dynamics (Verstrate et al. 2009, Krauter et 
al. 2011). While establishing rigorous robustness results requires further study, the ac- 
tual answer is expected to depend on the physical iniplenientation and its characteristic 
time scales. Lastly, the estimation of the speed of convergence still present numerous 
challenges, most importantly its optimization and a characterization of its scaling with 
the number of subsystems involved. 

Acknowledgements 

F.T. acknowledges support by the QUINTET and QFuture projects of the University 
of Padova. 

References 

[Alicki & Lcndi (1987)] Alicki, R. & Lendi, K. 1987 Quantum Dynamical Semigroups 

and Applications. Springer- Verlag, Berlin. 

[Barreiro et al. (2011)] Barreiro, J. T., MuUer, M., Schindler, P., Nigg, D., Monz, T., 
Chwalla, M., Honnrich, M., Roos, C. P., ZoUer, P. & Blatt, R. 2011 An open-system 
quantum simulator with trapped ions. Nature 470, 486-491. 

[Beige et al. (2000)] Beige, A., Braun, D., Tregenna, B. & Knight, P. L. 2000 Quantum 
computing using dissipation to remain in a decoherence-free subspace. Phys. Rev. 
Lett. 85, 1762-1765. 

[Bolognani & Ticozzi (2010)] Bolognani, S. & Ticozzi, F. 2010 Engineering stable 
discrete-time quantum dynamics via a canonical QR decomposition. IEEE Trans. 
Aut. Contr. 55, 2721-2734. 

[Combes et al. (2008)] Combes, J., Wiseman, H. M. & Jacobs, K. 2008 Rapid measure- 
ment of quantum systems using feedback control. Phys. Rev. Lett. 100, 160503:1-4. 

[Combes et al. (2010)] Combes, J., Wiseman, H. M. & Scott, A. J. 2010 Replacing 
quantum feedback with open-loop control and quantum filtering Phys. Rev. A 81, 
020301:1-4. 

[Diehl et al. (2008)] Diehl, S., MicheU, A., Kantian, A., Kraus, B., Biichler, H. P. & 
Zoller, P. 2008 Quantum states and phases in driven open quantum systems with 
cold atoms. Nature Phys. 4, 878-883. 

[Gorini et al. (1976)] Gorini, V., Kossakowski, A. & Sudarshan, E. G. C. 1976 Com- 
pletely positive dynamical semigroups of n-level systems. J. Math. Phys. 17, 821- 
825. 

[Hein et al. (2004)] Hein, M., Eisert, J. & Briegel, H. J. 2004. Multiparty entanglement 
in graph states Phys. Rev. A 69, 062311:1-20. 

[Khidjasteh et al. (2011)] Khodjasteh, K., Dobrovitski, V. V. & Viola, L. 2011 Pointer 

states via engineered dissipation. Phys. Rev. A 84, 022336:1-21. 

[Kraus et al. (2008)] Kraus, B., Diehl, S., Micheh, A., Kantian, A., Biichler, H. P. 
& Zoller, P. 2008 Preparation of entangled states by dissipative quantum Markov 
processes. Phys. Rev. A 78, 042307:1-9. 



10 



[Krauter et al. (2011)] Krauter, H., Muschik, C. A., Wasilewski, W., Petersen, J. M., 
Cirac, J. I. & Polzik, E. S 2011 Entanglement generated by dissipation and steady 
state entanglement of two macroscopic objects. Phys. Rev. Lett. 107, 080503:1-4. 

[Lindblad (1976)] Lindblad, G 1976 On the generators of quantum dynamical semi- 
groups. Commun. Math. Phys. 48, 119-130. 

[Lloyd & Viola (2001)] Lloyd, S. & Viola, L. 2001 Engineering quantum dynamics. 
Phys. Rev. A 65, 010101:1-4. 

[Perez-Garcia et al. (2007)] Perez-Garcia, D., Verstraete, P., Wolf, M. M. & Cirac, J. 
L 2007 Matrix product state representations. Quantum Inf. Comput. 7, 401-430. 

[Perez-Garcia et al. (2008)] Perez-Garcia, D., Verstraete, P., Wolf, M. M. & Cirac, J. 
L 2008 PEPS as unique ground states of local Hamiltonians, Quant. Inf. Comp. 8, 
0650-0663. 

[Pastawski et al. (2011)] Pastawski, F., Clemente, L. & Cirac, J. L 2011 Quantum 
memories based on engineered dissipation. Phys. Rev. A 83, 012304:1-12. 

[Poyatos et al. (1996)] Poyatos, J. P., Cirac, J. I. & ZoUer, P. 1996 Quantum reservoir 
engineering with laser cooled trapped ions. Phys. Rev. Lett. 77, 4728-4732. 

[Schirmer & Wang (2010)] Schirmer, S. G. & Wang, X. 2010 Stabilizing open quantum 
systems by markovian reservoir engineering. Phys. Rev. A 81, 062306:1-14. 

[Ticozzi & Viola (2008)] Ticozzi F. & Viola, L. 2008 Quantum Markovian subsystems: 
Invariance, attractivity and control. IEEE Trans. Autom. Contr. 53, 2048-2063. 

[Ticozzi & Viola (2009)] Ticozzi F. & Viola, L. 2009 Analysis and synthesis of attrac- 
tive quantum Markovian dynamics. Automatica 45, 2002-2009. 

[Ticozzi et al. (2010)] Ticozzi, F., Schirmer, S. G. & Wang, X. 2010 Stabilizing quan- 
tum states by constructive design of open quantum dynamics. IEEE Trans. Autom. 
Contr. 55, 2901-2905. 

[Ticozzi et al. (2011)] Ticozzi, P., Lucchese, R., Cappellaro, P. & Viola, L. 2011 Hamil- 
tonian control of quantum dynamical semigroups: Stabilization and convergence 
speed. IEEE Trans. Autom. Contr., in press. Preprint: I arXiv:llOl 24521 /1. 

[Verstraete et al. (2006)] Verstraete, P., Wolf, M. M., Perez-Garcia, D. & Cirac, J. I. 
2006 Criticality, the area law, and the computational power of projected entangled 
pair states, Phys. Rev. Lett. 96, 220601:1-4. 

[Verstraete et al. (2009)] Verstraete P., Wolf M. M. & Cirac, J. L 2009 Quantum com- 
putation and quantum-state engineering driven by dissipation. Nature Phys. 5, 633- 
636. 

[Wiseman & Milburn (2009)] Wiseman H. M. & Milburn G. J. 2009 Quantum Mea- 
surement and Control. Cambridge University Press, Cambridge. 



11 



